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Abstract 

By using the scaling method and the Thomas-Fermi and Extended Thomas- 
Fermi approaches to Relativistic Mean Field Theory the surface contribution 
to the leptodermous expansion of the finite nuclei incompressibility Ka has 
been self-consistently computed. The validity of the simplest expansion, which 
contains volume, volume-symmetry, surface and Coulomb terms, is examined 
by comparing it with self-consistent results of Ka for some currently used 
non-linear a — u parameter sets. A numerical estimate of higher-order con- 
tributions to the leptodermous expansion, namely the curvature and surface- 
symmetry terms, is made. 
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The curvature of the nuclear matter equation of state (EOS), i.e., the nuclear matter 
incompressibility is a key quantity in nuclear physics because it is related to many 
properties of nuclei (such as radii, masses and giant resonances), heavy-ion collisions, neutron 
stars and supernova collapses. One important source of information on is provided by the 
study of the isoscalar giant monopole resonance (GMR) (breathing mode) in finite nuclei. 
In the non-relativistic frame, theoretical microscopic calculations based on the random- 
phase approximation (RPA) fl[] and approximations to it such as the scaling method @-@] 
or constrained calculations J3|-^] using Skyrme |§ and Gogny || effective forces lead to a 
nuclear matter incompressibility coefficient of 215 ± 15 MeV |||7]]. A similar analysis 
carried out within the relativistic mean-field (RMF) theory with non-linear a — u effective 
Lagrangians gives a value of slightly higher: 250 — 270 MeV ||. 

The nuclear matter incompressibilty is not a directly measurable quantity, what is 
measured is, actually, the energy Em of the GMR of finite nuclei. It is convenient to write 
this energy in terms of the incompressibility for a finite nucleus of mass number A as: 



= V 1STT*P (1) 

where < r 2 > is the rms matter radius and M the nucleon mass. The finite nucleus in- 
compressibilty Ka can be parametrized by means of a leptodermous expansion [[| which is 
similar to the liquid drop mass formula: 

K A = K OD + KsfA" 1 ^ + K VS I 2 + K coul Z 2 A~^ 3 + , (2) 

where / = (N — Z) /A is the neutron excess. Eq. (H) suggests that it is possible to fit 
the coefficients of the expansion to the experimental data in a model independent way. 
Although some effort along these lines has been made in the past |J, the fact that a fit of 
the parameters of Eq. (0) to experimental data does not lead to a unique determination 
of the parameters is well established P, |l0| , |TT|| . Rather, the nuclear matter incompressibility 
has to be determined from effective forces which reproduce, in a microscopic calculation, the 
experimental values of the GMR excitation energy in heavy nuclei || . 



It is also possible to fit Ka calculated microscopically within the scaling model for a 
given effective interaction to the leptodermous expansion Eq. (|2|). This has been done, for 
example, in the non-relativistic frame using Skyrme forces |12[. In this case the coefficients 
entering Eq. (Q) can be expressed through infinite and semi-infinite nuclear matter prop- 
erties calculated with the Hartree-Fock approximation for each considered interaction. In 
particular, the volume-symmetry (K vs ) and Coulomb (K cou i) coefficients depend on some 
parameters of the liquid droplet model fl3| computed only using nuclear matter properties 



0. The surface coefficient K s f, also derived in 0, can be written as |L4 

54 



(3) 



The surface tension a is calculated in symmetric semi-infinite nuclear matter and is defined 



as: 



/+oo 
{H{p) - e O0 (p c )p}dz, (4) 
-co 

where p is the density profile whose central value is given by p c = p{— oo), Ti is the energy 
density and e ra is the energy per particle in nuclear matter at density p c . In Eq. (^) 
dots indicate the derivatives with respect to the central density and all the quantities are 
evaluated at a central density equal to the nuclear matter saturation density po, which is 
related to the radius constant r through 47rrQp /3 = 1. 

The key quantity entering Eq. ([3D is a which is the second derivative of cr(p c ) with 
respect to p c calculated at p c = p . The determination of a also requires knowledge of how 
the density profile p is modified during compression ||15|| . In the study of the breathing mode 



a scaling transformation of the densities is assumed. Actually, the coefficients entering the 
parametrization (0) can be derived under this hypothesis 0. The scaling transformation 
means that the density changes according to the transformation r —>■ Xr and consequently 

p x {r) = X 3 p(Xr). (5) 

Thus, in the scaling approach: 
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(6) 



To obtain the surface incompressibility coefficient K s f for a given effective interaction, it is 
necessary, first of all, to calculate the scaled surface tension o\ by replacing the densities 
by the scaled densities given by Eq. in Eq. ([|). In the non-relativistic frame this can 
be easily done within the Hartree-Fock scheme using zero-range Skyrme forces and a simple 



analytical expression for a\ is obtained |f2| . |i4 



The self-consistent calculation of K s f within the RMF approximation using the a — uj 
model is more involved due to the problem of the change in the meson fields induced by the 



scaled nuclear densities | 16|| . To our knowledge, only approximate calculations of K s j have 
been developed in the past for the relativistic model. This is the case of the Relativistic 
Thomas- Fermi (RTF) calculations of Refs. [fT6| , [r^] where a local density approximation of 
the meson fields was used. Another approach is related with the study of nuclei under 
an external pressure. Starting from a schematic energy density functional and adding a 
density-dependent constraint which simulates the pressure, analytical expressions for the 
surface tension a as a function of the bulk density p c can be derived for a wide class of 
compression modes, in particular, for the scaling mode This way one obtains the 

following formula for a in the scaling mode: 



Hpo) 



19 K^a 
81 po : 



(7) 



where a is the surface diffuseness parameter of a symmetric Fermi density. This pocket 
formula has been employed to estimate K s f in the RMF model for several non-linear o — uj 
parameter sets |L9| . A symmetric Fermi function that reproduces in the best way the density 
profile obtained from a Hartree calculation of semi-infinite nuclear matter has been used in 
ref . [^9| to determine the a parameter of Eq. (^) . 

Very recently, the scaling method applied to the RMF theory in the RTF and Relativistic 
Extended Thomas-Fermi (RETF) approaches has been used to self-consistently obtain the 
excitation energy of the GMR of finite nuclei |p0|,2l|. Our aim in the present paper is, firstly 



to obtain the surface coefficient K s f self-consistently in the RTF and RETF approaches 



developed in [pO| , pTf for some linear and non-linear a — lu parameter sets. On the other hand, 
we want to check whether the leptodermous expansion of the finite nucleus incompressibility 
Eq. (0) can reproduce the corresponding fully self-consistent value obtained in the RETF 
approach jn|] with some selected non-linear a — u parameter sets. 

The key point of our semi classical approach is that the local Fermi momentum kp and 
the effective mass m* scale as |j(J ,2j]]: 



kF\ = Xk F (Xr), rn*x( r ) = Am*(Ar), 



(8) 



where m* is still a function of A. With the help of Eq. (|§|), the nuclear part of the energy 
and the scalar density including h 2 corrections, which are functionals of k F and m*, scale 



as: 



S x (r) = A 4 £(Ar), p sX (r) = X d p a (Xr). 



(9) 



Again £ and p s are functions of the collective coordinate A because of their dependence on 
in*. Thus the scaled surface tension can be written as [p0|-|22| : 

.i d{Xz) 
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where all densities and fields depend on the variable Xz. With the help of the Klein-Gordon 
equations for the scaled vector and scalar fields derived from (|T0|), the scaled surface tension 
can be recast as: 



^ = I d(Xz) U 3 £x + \x 2 g v Vx P + ^g^xPs - ^ - ^ - A 2 eoo (A 3 P o)p 



Using the explicit RTF or RETF expressions for the nuclear part of the energy and for the 
scalar density p0|-p3| together with the Klein-Gordon equations for V\, <px, ^ and ^ 



derived from Eq. (|T0|), after some algebra the first and second derivatives of the scaled 



surface tension a\ with respect to A at A = 1 read (see Refs. p3| , |2lf for more details): 
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The first derivative of eoo(A 3 p ) at A = 1 is just three times the pressure calculated at satu- 



ration density and thus it vanishes, while the second derivative gives K^p pl| , p6| . On the 
other hand, since in the self-consistent RTF and RETF calculations the inputs for computing 
Eqs. (|l^)-(0) are quantities obtained from the solution of the variational equations associ- 
ated with the surface tension ( PH| ) at A = 1, the so-called "sigma dot" theorem is rigorously 
fulfilled [[TjJ . The method therefore allows a and consequently K s f to be computed on top of 
a self-consistent RTF or RETF calculation of the surface tension in symmetric semi-infinite 
nuclear matter. This is similar to what happens in the non-relativistic frame with Skyrme 



forces |14|, although in the relativistic case additional Klein-Gordon equations for and 



-r^ at A = 1 have to be solved. 



Now we shall discuss the results obtained from the self-consistent RTF and RETF meth- 
ods in the scaling approximation. Table 1 collects K^, a and K s f for the non-linear NL-Z2 
2§, NL1 @,NL3 |0|, NL-RA1 |3lf, NL-SH |2f and NL2 |3f and the linear HS f3| and 
LI parameter sets. One observes that in both the RTF and RETF calculations a and 
K s f decrease (become more negative) with increasing bulk incompressibility K^. The RTF 
and RETF values of a and K s f for a given parameter set are, in general, rather different 
from one another, which means that the precise value of these quantities is model dependent. 
This is known to happen also with other quantities related with the nuclear surface. For 
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example, such is the case of the surface energy coefficient of the leptodermous expansion of 
the binding energy of a nucleus, which is calculated as Airrfo. The quality of the RTF and 
RETF approximations for semi-infinite nuclear matter and finite nuclei with respect to the 
RMF Hartree approach, and its dependence on the effective interaction, was investigated 
in Refs. P2],E3,PS[ by analyzing the results obtained with many different parameter sets. It 



was noticed that the RTF results fluctuate around the corresponding Hartree results. Due 
to this fact there exist parametrizations for which the RTF approximation agrees by chance 
with the Hartree result. The behaviour of RETF results in comparison with the Hartree 
solutions was found to be less dependent on the parameters of the force than in the RTF 
case, and it turned out that on the average the RETF results are in better agreement with 



the Hartree ones 22,24.25 



The first contribution to K s f in Eq. (|3]) comes from the surface tension, let us call it 
K°f. The deviation found in the value of the surface tension from RTF calculations with 
respect to the corresponding RMF Hartree calculations is strongly correlated with the value 
of the effective mass in nuclear matter m^/m J2^,^5|]. For small values of mj^/m the RTF 
surface tension is larger than the Hartree one, it practically agrees with the Hartree result 
for rrf^/m ~ 0.65, and it becomes smaller than the Hartree result for larger m^/m. On the 
other hand, the RETF result for the surface tension is consistently lower than the Hartree 
result and much less dependent on the specific value of m^/m. (A similar situation is found 
for the total energy of finite nuclei p^ , |24]p5| .) These trends, of course, are also reflected 
in K a sf . For example, for NL1 (m^/m = 0.57) we have K a sf = 402.6, 377 and 429.3 MeV 
in the Hartree, RETF and RTF approaches, respectively. For NL2 (m^/m = 0.67) it is 
Kf f = 479.6, 439.1 and 465.7 MeV in the Hartree, RETF and RTF calculation, respectively. 

The second contribution to K s f in Eq. (|3]) is due to the second derivative of the surface 
tension. The results for a in the RTF approach decrease with faster than in the RETF 
calculation. At small values of the RTF value of a is less negative than that computed in 
the RETF approach, while the opposite happens for higher values of K^. Both approaches 
predict the same value of d for an incompressibility around that of NL1 (211 MeV). A 



similar behaviour is displayed by the self-consistent values of K s f, although the crossing 
point between the RTF and RETF predictions is shifted to a larger value of (around 
280 MeV) due to the fact that the contribution proportional to o (K°j) is larger in the RTF 
approach than in the RETF approach for the parameter sets considered here. 

The a values obtained from the pocket formula Eq. (0) using the surface diffuseness 
of the RTF or RETF semi- infinite nuclear matter density profiles also decrease with K^, 
though the estimate provided by Eq. (0) does not reproduce in general the self-consistent 
values very precisely. The approximate a is always smaller when calculated from the RTF 
approach than from the RETF approach for the parametrizations of Table 1. Using Eq. (|7|) 
to estimate the value of K s f in first approximation, one finds that the RETF result is closer 
to the Hartree value than the RTF result for the sets NL-Z2, NL1, NL3 and NL-RAl. For 
NL2 and LI it is the RTF estimate which lies closer to the Hartree estimate. For NL-SH 
and HS the approximate Hartree prediction lies roughly in between of the approximate RTF 
and RETF values. To the extent that Eq. (0) is applicable, it provides a hint of where the 
unknown exact Hartree value for K s f should lie with respect to the self-consistent RTF and 
RETF results presented in Table 1. 

Another different approach to computing K s f was proposed in Refs. |TI||T7|] . It is based 
on the scaling method together with a local density approximation for the meson fields 
within the RTF approach. In Ref. |17j a K s f of approximately —1000 MeV was reported for 
a linear set with = 545 MeV (similar to the HS set). This result is in good agreement 
with that of our self-consistent RTF calculation for the HS set, though it is clearly larger 
in absolute value than the RETF result for K s f obtained with the HS set. On the other 
hand, our self-consistent RTF and RETF surface incompressibilities differ considerably from 
the estimate of Ref. where approximate K s f values of —333.1 and —610.1 MeV were 



reported for the NL1 and NL-SH parametrizations calculated with the method used in [jTo 

It should also be pointed out that in our self-consistent semiclassical calculations we 
find that the ratio between the surface and bulk incompressibilities increases with (in 



agreement with the results of [0]). In the RETF case this ratio is close to one, as happens 
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for the non-relativistic Skyrme forces ||, provided that the bulk incompressibility of 
the interaction is not excessively high. In the RTF case the ratio between the surface and 
bulk incompressibilities increases much faster with Koo than in the RETF calculations, and 
it considerably differs from unity for parametrizations with either a very low or a very high 
bulk incompressibility. In Figure 1 we plot —K s f as a function of K^ for the parameter 
sets considered in Table 1. As in the non-relativistic case K s f varies roughly linearly 
with Koo. A linear fit of all the points gives —K s f = lAlK^ — 84 in the RETF model 
and —K s f = 2.19/Coo — 295 in the RTF model. If only the non-linear parametrizations are 
included in the fit one obtains —K s f = 1.35Koo — 54 and —K s f = 1.96/Coo — 238 in the 
RETF and RTF cases, respectively. 

The surface incompressibility coefficient is both large and negative, thus its contribution 
considerably reduces the finite nucleus incompressibility Ka with respect to the nuclear 
matter limit K^. This result, although obtained in the scaling model, illustrates the physical 
effect that the compression of the surface provides a considerable reduction of Ka, which is 
also found in more fundamental RPA calculations 0. In Ref. [^T|] we have self-consistently 
computed the finite nucleus incompressibility Ka using the RETF approach and the scaling 
method which we have employed in the present work to obtain K s f. Thus we can now 
precisely check the ability of the leptodermous expansion Eq. @ in reproducing the full 
calculation of Ka carried out in Ref. 0] for various finite nuclei. 

The coefficients K vs and K coui entering Eq. (H) are computed using nuclear matter 
properties only. Explicit expressions for these coefficients in the non-linear a — uj model 
are reported in fL8| . In our analysis we will use the NL1, NL3 and NL-SH parameter 
sets for which the numerical values of these coefficients are given in [[HJ. The surface 



incompressibility coefficient is the self-consistent value taken from Table 1. Table 2 collects 



Ka obtained from the full self-consistent RETF calculation |2]J as well as the value K(A, I) 
given by Eq. <g) for 40 Ca, 48 Ca, 56 Ni, 90 Zr, 116 Sn, 132 Sn, 144 Sm and 208 Pb. From this Table it 
can be seen that the leptodermous expansion with the terms given in Eq. fails to describe 
small nuclei and also very asymmetric nuclei such as 132 Sn or 208 Pb. In addition, some words 
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of caution should be said about the Coulomb term in Eq. (0). In the self-consistent scaling 
calculation of the finite nucleus incompressibility, the Coulomb energy does not participate 
directly if the scaling Eq. (|5|) for the density is assumed to apply [|3|,f2T| . Thus, the Coulomb 
term in Eq. (fj) should be related to the change in Ka when the Coulomb interaction is 
switched off in the self-consistent calculation. The Coulomb term in Eq. (|2[) overestimates 
this change by approximately 6 MeV for NL1, 3 MeV for NL3 and 1 MeV for NL-SH. 

Now we would like to analyze whether the addition of some higher order terms in the 
leptodermous expansion Eq. (fj) improves the agreement with the Ka results calculated self- 
consistently. In particular, we will focus our attention on the curvature K CV A^ 1 ^ and the 
surface-symmetry K SS I 2 A~ 1 ^ terms. Although these terms could be derived by enlarging 
the leptodermous expansion of Blaizot ||, as has been done in the non-relativistic case 
E| , it becomes more complicated in the relativistic case. Thus, for a fast estimate of the 



curvature and surface-symmetry terms, we perform a numerical fit. To do this, we follow 



the same strategy as in Ref . (12| . First we consider symmetric nuclei with the Coulomb force 
switched off. In this case the leptodermous expansion Eq.(^j) (adding the curvature term) 
reduces to: 

K A = K 00 + K sf A~ 1 ^ + K CV A~ 2 '\ (14) 

In Figure 2 we plot [Ka — K^A 1 ^ 3 versus A~ 1 ^ for the three parameters sets used in this 
analysis. Here K^ is the nuclear matter incompressibility given in Table 1 and Ka are the 
self-consistent incompressibilities calculated for A ranging from 300 up to 300000. In the 
linear fit of these curves the y— axis intercept gives K s f of the corresponding force, while 
the slope gives K cv . The surface terms obtained in this way are —246.1, —328.4 and —435.8 
MeV for the NL1, NL3 and NL-SH parameter sets, which are very close to the corresponding 
self-consistent values (see Table 1). The estimates of the curvature term in the leptodermous 
expansion of the finite nucleus incompressibility obtained with NL1, NL3 and NL-SH are 
-317.2, -229.8 and -185.6 MeV respectively. 

To obtain the surface-symmetry contribution, we have found it convenient to parametrize 
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the difference between the self-consistent incompressibilities Ka of a given nucleus with 
neutron excess I and the corresponding symmetric nucleus as: 

K A ,i - K A ,i=o = K VS I 2 + K^A- 1 '^ (15) 

where again uncharged nuclei have been considered. For each parameter set and accord- 
ing to (|15|), if [Ka,i — Ka,i=q]I~ 2 is plotted versus A^ 1 ^ a unique curve should be found 
which is independent of the value of I. However, one obtains a family of almost parallel 
lines whose slope is K ss . The splitting of these lines gives us information on the higher 
order symmetry contributions missed in the parametrization ([15]). Thus we will estimate 
the surface-symmetry term from a linear fit of the curve corresponding to / = 0.1, which 
roughly corresponds to an average asymmetry along the periodic table. This curve is plotted 
in Figure 3 for A ranging from 250 to 200000 for each considered parameter set. The corre- 
sponding y-axis intercepts agree very well with the K vs values calculated in nuclear matter 
(-676.1, -698.9 and -794.5 MeV for NL1, NL3 and NL-SH respectively @). Our estimate 



of the surface-symmetry contribution to Ka corresponds to the slopes of these linear fits, 
which are 1951.4, 1754.0 and 1716.5 MeV for NL1, NL3 and NL-SH respectively. 

Table 3 collects the self-consistent finite nuclei incompressibility Ka (without Coulomb) 
compared with the macroscopic parametrizations K (A, I) (Eq. (§)) and K*(A,I) which 
contains the curvature and surface-symmetry contributions obtained from the previously 
discussed fits. Again, the self-consistent incompressibilities corresponding to the lightest 
nuclei and the very asymmetric nuclei are not well reproduced by the simplest expansion 
Eq. (0). If the curvature and surface-symmetry corrections are included, the improved 
macroscopic formula K*(A, I) reproduces the self-consistent incompressibilities with an er- 
ror, on average, smaller than 1.2%, 0.9% and 0.3% for the NL1, NL3 and NL-SH parameter 
sets. In order to gain some insight into the accuracy of our estimate of the curvature and 
surface-symmetry contributions, we fit the self-consistent results for the finite nuclei con- 
sidered in Table 3 to a leptodermous expansion including curvature and surface-symmetry 
terms. The volume, surface and volume symmetry coefficients are taken from self-consistent 
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infinite and semi-infinite nuclear matter calculations. The results of this calculation show 
that the difference of the curvature contribution obtained from the fit in the asymptotic 
region and from finite nuclei is always less than 10%, whereas the difference in the surface- 
symmetry contribution lies below 3%. 

We have applied the scaling method in the Thomas-Fermi and Extended Thomas-Fermi 
approximations to the Relativistic Mean Field Theory to self-consistently calculate the sur- 
face coefficient K s f of the leptodermous expansion of the finite nucleus incompressibility 
derived within the Blaizot model. The ratio between the surface and bulk incompressibilities 
obtained in our semiclassical calculation increases with the nuclear matter incompressibility, 
more strongly in the RTF than in the RETF case. In the RETF calculations this ratio is 
close to one, as in the case of non-relativistic Skyrme forces, for the non-linear parameter 
sets which have a nuclear matter incompressibility not larger than roughly 300 MeV. 

For the analyzed a — uo parameter sets, the leptodermous expansion Eq. (|2|) is not able 
to reproduce very precisely the finite nuclei incompressibilities obtained self-consistently. 
In particular the macroscopic contribution of the Coulomb force can differ from the self- 
consistent contribution up to 6 MeV. We have numerically estimated higher order contri- 
butions to the leptodermous expansion, namely curvature and surface-symmetry terms, in 
the asymptotic region (i.e., for very large uncharged systems). We have found that the fi- 
nite nuclei incompressibilities are reasonably well reproduced by an extended leptodermous 
expansion which includes curvature and surface-symmetry contributions. 
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FIGURE CAPTIONS 



Figure 1. Surface incompressibility coefficient versus the nuclear matter incompressibility 
modulus for the parameter sets of Table 1. 

Figure 2. (K A — K^A 1 / 3 versus A^ 1 ^ computed for several uncharged and symmetric 
nuclei from A = 250 to A = 300000 for the NL1, NL3 and NL-SH parameter sets. 

Figure 3. (Ka,i — Ka,i=o) 1 1 2 versus A^ 1 ^ 3 for several uncharged nuclei from A = 250 to 
A = 200000 with a neutron excess 0.10 for the NL1, NL3 and NL-SH parameter sets. 
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TABLES 

TABLE I. Values of a (in MeVfm 4 ) and K sf (in MeV) calculated with the RTF and RETF 
approaches and the scaling method for several parameter sets. The nuclear matter incompressibility 
modulus Koo (in MeV) and the —Kgf/Koo ratio are also listed. 









RTF 






RETF 




ij 


Ksf 


—Ksf /Koo 


o 


Ksf 


-Ksf/K^ 


NL-Z2 


172.2 


— 113.9 


-85.2 


0.49 


— 131.2 


-182.5 


1.06 


NL1 


211.1 


— 170.3 


-170.6 


0.81 


— 171 .8 


-225.4 


1.07 


NL3 


271.5 


-224.2 


-310.4 


1.14 


-209.3 


-313.7 


1.16 


NL-RA1 


285.3 


-235.5 


-335.4 


1.18 


-216.6 


-326.7 


1.15 


NL-SH 


355.0 


-292.7 


-469.8 


1.32 


-258.2 


-429.6 


1.21 


NL2 


399.2 


-295.9 


-521.0 


1.31 


-279.0 


-482.8 


1.21 


HS 


546.8 


-521.5 


-996.7 


1.82 


-424.9 


-804.2 


1.47 


LI 


625.6 


-422.6 


-1024.6 


1.64 


-320.6 


-787.1 


1.26 
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TABLE II. Finite nuclei incompressibilities (in MeV) calculated with the self-consistent RETF 
approach (Ka) and with the leptodermous expansion Eq. (g) (K(A, I)). Results are presented for 
the NL1, NL3 and NL-SH parameter sets. 





NL1 




NL3 




NL-SH 






K A 


K(A,I) 


K A 


K(A,I) 


K A 


K(A,I) 


40 Ca 


108.2 


128.1 


145.3 


161.0 


196.8 


208.6 


48 Ca 


111.1 


116.9 


147.4 


151.0 


198.3 


198.4 


56 Ni 


115.0 


130.8 


153.2 


166.0 


207.1 


216.7 


90 Zr 


122.5 


129.3 


161.6 


167.3 


217.5 


221.1 


116 Sn 


124.3 


126.3 


163.4 


165.4 


219.8 


220.4 


132 Sn 


121.3 


105.4 


157.6 


144.9 


210.9 


197.5 


144 Sm 


125.4 


125.3 


164.5 


165.3 


221.6 


221.5 


208p b 


124.1 


111.1 


161.1 


152.1 


216.7 


208.1 
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TABLE III. Finite nuclei incompressibilities (in MeV) for several uncharged nuclei calcu- 
lated self-consistently using the RETF approach (Ka), with the leptodermous expansion Eq. (|2|) 
(K(A, I)) and including the curvature and surface-symmetry contributions (K*(A, I)). Results are 
presented for the NLl, NL3 and NL-SH parameter sets. 







NLl 






NL3 






NL-SH 






K A 


K(A, I) 


K*(A,I) 


K A 


K(A,I) 


K*(A, I) 


K A 


K(A, I) 


K*(A,I) 


40 Ca 


118.6 


145.2 


118.1 


160.1 


179.8 


160.2 


213.4 


229.4 


213.5 


48 Ca 


119.6 


130.3 


121.2 


159.7 


165.8 


161.8 


215.1 


214.7 


213.8 


56 Ni 


129.3 


152.2 


130.5 


172.7 


189.6 


173.9 


230.5 


242.7 


230.0 


90 Zr 


139.6 


152.5 


142.0 


184.2 


192.9 


186.3 


244.3 


249.3 


244.8 


116 Sn 


144.0 


152.0 


146.3 


189.0 


193.9 


191.1 


250.1 


251.8 


250.7 


132 Sn 


137.2 


127.1 


137.4 


179.3 


168.9 


180.2 


236.8 


223.9 


236.58 


144 Sm 


148.5 


155.0 


150.7 


194.3 


198.2 


196.3 


256.7 


257.7 


257.2 


208p b 


148.4 


142.8 


148.5 


193.4 


187.3 


194.0 


255.0 


246.9 


254.6 
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